The limit of strong ion coupling due to electron shielding 
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We show that strong couphng between ions in an ultracold neutral plasma is limited by electron 
screening. While electron screening reduces the quasi-equilibrium ion temperature, it also reduces 
the ion-ion electrical potential energy. The net result is that the ratio of nearest-neighbor potential 
energy to kinetic energy in quasi-equilibrium is constant and limited to approximately 2 unless 
the electrons are heated by some external source. We support these conclusions by reporting new 
measurements of the ion velocity distribution in an ultracold neutral calcium plasma. These results 
match previously reported simulations of Yukawa systems. Theoretical considerations are used to 
determine the screened nearest-neighbor potential energy in the plasma. 
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I. INTRODUCTION 

Ultracold neutral plasmas are strongly coupled 
Coulomb systems. Strong coupling occurs vifhen the 
nearest-neighbor (Coulomb) potential energy is greater 
than the average kinetic energy per particle. The degree 
of strong coupling between ions is traditionally charac- 
terized by the parameter 
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where Z is the ion charge state, q is the fundamental 
charge, a^g = (3/4?™)^/^ is the Wigner-Seitz radius or 
the average distance between ions in the plasma, n is 
the plasma density, and T is the temperature. These 
plasmas are typically generated by photo-ionizing laser- 
cooled gases [1] or ions in an ultrasonic jet [2]. They are 
diagnosed using three-body recombination [3-8], ther- 
malization rates [9-12], electron evaporation or rf absorp- 
tion [3, 13-18], charged particle imaging and detection 
[2, 19], and optical fluorescence [12, 20, 21] and absorp- 
tion [22-24]. Theoretical calculations and simulations 
[25-32] give great insights into the properties of these 
plasmas. 

Achieving a higher value of the strong coupling pa- 
rameter is a high priority in this field. One of the pri- 
mary objectives of ultracold plasma studies is under- 
standing how strong coupling changes the description of 
basic plasma or atomic processes. As the plasma be- 
comes more strongly coupled, the fundamental assump- 
tions that make fluid approximations valid break down. 
Collisional processes that are well-understood in terms 
of interactions between discreet atoms or ions in the 
plasma, such as three-body recombination, should take 
on a many-body nature. While evidence for this is seen 
for dense Rydberg gases [33], only limited evidence for 



departures from traditional plasma physics has been re- 
ported in the literature (see, for example [5, 10]). 

At early times in the plasma evolution the ion temper- 
ature is dominated by disorder- induced heating (DIH). 
Although ions are formed with essentially zero kinetic 
energy in ultracold neutral plasmas, they have an excess 
of electrical potential energy due to the random spatial 
distribution of nearest neighbors in the plasma. In the 
absence of correlation and shielding effects in the final 
state, the initially T ~ mK ions heat up to the correla- 
tion temperature 
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on the time scale of the inverse ion plasma frequency, 
l/cjj (ng^/TOjeg)""'^/^, where to, is the ion mass. In 
the absence of electron shielding, the ion temperature is 
determined by the density alone. 

In neutral plasmas, electrons form a screening back- 
ground for the ions. This changes the ion dynamics. If 
the electron temperature is not too low, the ion interac- 
tion can be modeled with a Yukawa potential 
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where Ad — {k^T^eQ/nq^Y/'^ is the Debye length and 
fcg is Boltzmann's constant. Electron screening can be 
parameterized using the inverse scaled screening length, 
K = a^g/Xj^. The parameter k exhibits a strong tem- 
perature dependence, k ^ Te and a weak density 
dependence k. ~ n^/^ . 

Previous studies have shown that electron screening 
can significantly reduce the ion temperature and slow 
the DIH time scale [10, 12]. At first sight this appears to 
increase T. However, because electron screening reduces 
the nearest-neighbor electrical potential energy, the net 
effect of electron shielding on the ratio of the actual 
nearest-neighbor potential energy to kinetic energy is not 
immediately clear. 
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In this paper we present new fluorescence measure- 
ments of an ultracold neutral Ca plasma. From our 
measurements we extract the time evolving rms width 
of the velocity distribution by fitting the data to a Voigt 
profile. Using an expansion model we find the electron 
and ion temperature as a function of time. By varying 
the initial electron temperature we generate plasmas with 
varying degrees of electron shielding. We show that we 
can generate plasmas with very cold ions by mitigating 
the effects of DIH through electron shielding. We com- 
pare our experimental results to two molecular dynam- 
ics simulations, which show good agreement with each 
other and with our data. We use theoretical considera- 
tions to extract the screened ion-ion potential energy in 
the plasma. We suggest that although electron screening 
reduces heating due to DIH, it also reduces the nearest- 
neighbor potential energy in such a way that the ratio of 
potential energy to kinetic energy is independent of the 
electron temperature. 



II. EXPERIMENT 

We cool and trap approximately 20 million "^^Ca atoms 
in a magneto-optical trap (MOT) at a temperature of 
a few mK. The spatial density profile is Gaussian and 
has the form n(r) = Uq exp(— r^/2cr^), with peak density 
rig < 2 X 10-^° cm~^ and a = 0.2 - 0.4 mm. We ionize 
these laser-cooled atoms using ns-duration laser pulses at 
wavelengths of 423 and 390 nm that drive the 4s^ ^Sq ^ 
4s4p and the 4s4p — >■ continuum transitions, 
respectively. 

The initial ion and electron temperatures, Tj(G) and 
Tg(G), are controlled in the experiment but change in 
time. The initial ion temperature in the plasma is ap- 
proximately equal to the few mK temperature of the neu- 
tral atoms in the MOT. It changes as ions undergo DIH 
and as the plasma expands. The plasma is not confined 
by the MOT, and it freely expands into the surrounding 
vacuum. The expansion velocity is typically determined 
by electron temperature and the ion mass. As long as the 
electron temperature is not too low, the asymptotic ex- 
pansion velocity is v = {k^Tjm^y/"^ [34]. The initial 
electron energy in the plasma is equal to the difference 
between the ionizing laser photon energy and the calcium 
ionization energy. In our experiment this typically ranges 
from = 2Ej3kB = 0.5 - 150 K. 

We probe the ion velocity distribution using laser- 
induced fluorescence. We use a low-power cw laser beam 
detuned from resonance at a wavelength of 397 nm, cor- 
responding to the Ca"*" 4s ~^ 4p '^P°/2 transition. 
The laser beam is coUimated to a diameter of 4 mm, at- 
tenuated to 2.5 mW, aligned to spatially overlap the cal- 
cium plasma, and retroreflected. The maximum probe 
laser beam intensity in the retroreflected conflguration 
is approximately equal to the saturation intensity, where 
the saturation intensity is Igat = 46 mW/cm^. Fluores- 
cence photons at this same wavelength are collected as 



a function of time after the plasma is generated using a 
1-GHz bandwidth photo-multiplier tube and digital os- 
cilloscope. We repeat this process for a range of probe 
laser beam offset frequencies from to about ±250 MHz. 



III. DATA ANALYSIS 

We extract the time evolving ion velocity ^^^^ by flt- 
ting the fluorescence data to a Voigt profile. The Voigt 
profile is a mathematical representation of the absorp- 
tion cross section per atom. It is the convolution of a 
Lorentzian and a Gaussian lineshape 
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with the Lorenztian and Gaussian profiles given by 
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respectively. In these equations, v is the detuning from 

resonance, 7 is half the natural linewidth of the atomic 
transition (the HWHM of the Lorentzian line shape) , and 
t'j.^g is the rms Gaussian width. The integral defined in 
Eq. (4) can be evaluated as 



Re[u;(^;)] 



2-KVr 



(7) 



The term in the numerator is the complex error func- 
tion, and it is given by w{z) = erfc(— 10), where z 
is {v -I- i7)/-\/2fj.jjjg, and erfc is the complementary er- 
ror function. In the analysis, the Lorenztian half width 
is equal to 11 MHz, half the natural linewidth of the 
397 nm transition. However it is power-broadened to 
13.2 MHz (HWHM). The Gaussian width v^^^ is ex- 
tracted as a fit parameter. It is converted to the rms 
width of the velocity distribution using the Doppler shift, 

«i,rms = (^B^j/™i)^^^ = •^^'rms^ allowing US to map out 
the width of the ion velocity distribution as a function of 
time. 

As the plasma evolves, the plasma ions are accelerated 
radially outwards. The acceleration is well approximated 
by the expression a ~ {k^Te/mi(j'^)r [25]. This directed 
expansion adds to the thermal motion of the ions. In our 
experiment we measure the total rms width of the veloc- 
ity distribution, which contains contributions from both 
the random thermal motion as well as the accelerated 
expansion. Fortunately, the time scale for expansion is 
slower than the time scale for DIH, and the contributions 
of each of these effects can be cleanly separated. 
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FIG. 1. The time evolving rms width of the ion velocity 
distribution for an ultracold calcium plasma at a density of 
approximately 0.9 x 10^" cm^'^ and electron temperature of 
Te = 50 K. The rms velocity width is found using a fit to a 
Voigt profile, where the Gaussian frequency width is extracted 
as a fit parameter and converted to the velocity width through 
the Doppler shift. The model described by Eq. (8) is plotted 
as the gray dashed line. 



FIG. 2. The equilibrium ion temperature after DIH (trian- 
gles) and the characteristic DIH time (circles) plotted as a 
function of the electron temperature extracted from the ex- 
pansion model. It is evident from this plot that electron 
screening reduces the ion temperature and extends the DIH 
time. The density for all these plasmas is approximately 
0.9 X 10^° cm"^ 



Using the found from the Voigt fitting we are 

able extract the ion temperature. The ion temperature 
can be related to the rms ion velocity as [21] 

y "''i I ' exp J 

where Texp, the characteristic expansion time, is given 
by Texp = Vm,a{0)ykB[T,{0)+Ti{0)] and the time 
evolving ion and electron temperatures are given by 
T^(t) = T'«(0)/(1 + t^/Tgxp) and the subscript a = i,e. 
In Fig. 1 we have plotted the v^ calculated from this 
model and the ^^^^ found by fitting a Voigt profile to 
the fluorescence data of one of our plasmas. Rearranging 
Eq. (8) we solve for the ion temperature Ti{t) 

^exp 

The initial ion temperature is assumed to be negligible 
compared to the initial electron temperature. On the 
200 - 1000 ns timescale the plasma has not expanded. 

When the electron temperature is not known, it is pos- 
sible to estimate its value using the expansion model 
in Eq.(8). This is useful for plasmas that evolve from 
cold Rydberg gases, for example, where the initial elec- 
tron temperature is not well defined. At these small 
initial electron temperatures, the electron temperature 
extracted from the expansion model have been shown to 
overestimates the electron temperature at late times [34]. 
However, it is likely that the model is reasonably accurate 
at early times before the plasma has expanded. 



IV. EXPERIMENTAL RESULTS 

As the initial electron temperature decreases, the 
plasma expansion rate and ion thermalization rate 
changes. In Fig. 2 we plot the average T^ = m^v'^^^/k-Q 
after the DIH process has completed. We also plot the 
characteristic DIH time for a range of electron temper- 
atures. This characteristic time is taken to be the time 
when the ion temperature oscillation due to DIH is at 
its maximum. This can be seen in Fig. 1 at approx- 
imately 110 ns, where the measured t^i^ms (black line) 
passes slightly above the model (gray dashed line). The 
electron and ion temperatures arc both extracted from 
the expansion model given by Eq. (8). Plasmas with 
smaller T,. have smaller values of . The electron shield- 
ing length deceases with decreasing T^, softening the ion- 
ion Coulomb interaction. As the ions move under the 
influence of the screened Coulomb force of the neighbor- 
ing ions, they acquire less kinetic energy compared to the 
unscreened case. 

As decreases, the time scale for DIH increases. This 
confirms observations in Refs. 10 and 12. The low- 
temperaturc electrons more effectively shield ions from 
their nearest neighbors. The Coulomb force is reduced, 
and the ions take longer to reach their "equilibrium" po- 
sitions. The data in Fig. 2 show that the DIH time is 
extended by as much as a factor of 2. 



V. COMPARISON WITH MODELS 

Molecular dynamics simulations of complex neutral 
plasmas were published in the 1990's by Farouki, Ham- 
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aguchi, and Dubin [35-38]. Those simulations showed 
that electron shielding and correlation effects reduce the 
average electrical potential energy of the plasma ions. 
Murillo showed that these simulations can be applied to 
ultracold neutral plasmas [39] . A discussion of how that 
is done is given below. 



A. Deriving the DIH ion temperature 

The energy density per particle of the ultracold plasma 
can be generically written as 
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where the siimmation indexes represent electrons e or 
ions i. The Coulomb potential is written as u^;, = 
(g„g^/47reQ)(l/r) and the radial distribution function be- 
tween species a and b is gai,{r). This expression assumes 
a uniform plasma density, n. While the potential energy 
terms in Eq. (10) are general, the kinetic energy is writ- 
ten in terms of the temperatures and T!- . 

In order to quantify the ion heating during the DIH 
process, we will examine Eq. (10) at two important in- 
stances in the plasma evolution. The first instance is just 
after the plasma is formed, after the electrons have ther- 
malized with each other but before the ion have moved 
(1/Wg ~ 1 ns). Compared to the other energy scales in 
the system at this moment, the ions have essentially zero 
kinetic energy and we will therefore set the to zero. 
We will call this instance the initial time. 

The second instance is after the ions have thermal- 
ized with each other but before the plasma has expanded 
(1/wj ~ 100 — 500 ns). The ions will have moved pri- 
marily due to the Coulomb force of their nearest neigh- 
bors. Collisional transfer of energy from the electrons to 
the ions is much slower than the ion-ion collision rate. 
Consequently the electrons and ions maintain separate 
temperatures. The kinetic energy gained by the ions to 
this point will have come from the electrical potential en- 
ergy of their screened neighboring ions. We will call this 
instance the final time. 

Conserving energy, £j = £p, allows us to write. 
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where £j and £p are the initial and final energy densi- 
ties and the explicit r-dependence of the potentials and 
distribution functions has been suppressed. 

We can solve this equation if we assume that the elec- 
tron temperature remains constant. In this case we can 



also ignore the term u^^g^^ on both sides of the equation 
because it remains unchanged from the initial to final 
time. We can also ignore the term u^^g^^ because the 
constant approximation doesn't change the coupling 
between the electrons and ions when the ions move. The 
dominant change occurs in the ion-ion interaction. 
With these approximations, Eq. (11) becomes 
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Note that the term g^^ in the initial state has been 
dropped because the initial state is completely disordered 
and [^jjj = 1. We can make a connection with the 
Yukawa-MD simulations of Refs. 35-38 by introducing 
the Yukawa potential. The final state ion-ion potential 
can be trivially written as 

<i\ ■ (13) 

Similarly, we can express the radial distribution function 
as 

9u{r) = [guir)-l] + l = h,,{r) + l, (14) 

where h^^{r) is the pair correlation function. Inserting 
these definitions into Eq. (12) and simplifying gives 
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where all quantities are evaluated in the final state. The 
second term on the RHS is small and can be neglected. 



At small r, the quantity u^. 



is small and at large r 



the pair correlation function h,i^{r) tends to zero. The last 
term on the RHS can be evaluated directly and is equal to 
-{Z/2)nTkBTjK^. The first term on the RHS has been 
tabulated using molecular dynamics (MD) simulations 
[35 38]. It is the potential energy of the Yukawa ions 
after the DIH process has completed. In order to compare 
directly with the MD simulations, we need to convert 
from energy density to energy per particle. This is done 
by multiplying by the volume and dividing by the number 
of ions. We find the final ion temperature to be 



l^sT, = -| / d'rulih,, + 1) - IJ^fc^T,. 



(16) 



While this expression could be simplified further, we will 
leave it in this form in order to more easily compare with 
the results of previously published MD simulations. 

The MD simulations by Hamaguchi tabulate the 
temperature-scaled "excess energy" per particle, u = 
U/NkT, which is written as [37] 
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We recognize the first term on the RHS of Eq. (17) as the 
integral in Eq. (16) divided by k^T^. The last term on 
the RHS of Eq. (17) explicitly accounts for the periodic 
boundary conditions, which we will neglect because we 
are considering an infinite-sized plasma. Equation (17) 
includes the energy of the Debye sheath, —k/2. To get 
the "true" potential energy per particle, we add this back 
in and multiply by k^T^, 

N — l N 
j=l k=j+l 

(18) 

Comparing this with Eq. (16) gives 
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This derivation is complementary to the one presented in 
[39]. It's presentation here provides information regard- 
ing the approximations and assumptions used in obtain- 
ing this result. 



B. Experimental and theoretical DIH ion 
temperature 

Using Eq. (19) and the MD results of Hamaguchi et al. 
[35-38] , we can predict the ion tcrnpciraturc after the; DIH 
process has completed. This determination requires an 
iterative process because the ion temperature appears in 
the RHS of Eq. (19) in the u/T term [23]. One begins by 
choosing an initial ion temperature, density, and electron 
K. From this the ion F can be calculated. The tables 
in Refs. 35-38 then give the value of u/T. This can be 
inserted into Eq. (19) which gives a new ion temperature. 
The process is repeated until the ion temperature and F 
converge to a self-consistent limit. The resulting F as a 
function of k is plotted in Fig. 3. As k increases, the 
Debye length Xjj becomes smaller. The electrons more 
effectively shield the neighboring ions from one another 
and the final DIH temperature decreases. The F vs. k 
plot is a plot of l/T, vs. 

A more recent MD simulation was published that cal- 
culated the influence of electron shielding on the ion 
DIH temperature [10]. This MD simulation of Yukawa- 
shielded calcium ions started with the ions at rest, ran- 
domly positioned in a cell. The ions were allowed to move 
in the field generated by all of the other shielded ions in 
the cell (with periodic boundary conditions). After sev- 
eral u)^^, the average ion kinetic energy was calculated 
and from this the ion temperature was determined. The 
result is givcm in Eq. (6) in Ref. 10. This is plotted 
as a dashed line in Fig. 3. The k domain of this calcu- 
lation is somewhat limited because the plasma becomes 
non-ideal when k > 1 in the singly-ionized plasma used 
in the study, and in those conditions for that plasma it 
is not clear that the Yukawa approximation is valid. The 
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FIG. 3. Theoretical and oxporimental plots of the coupling 
parameter F as a function of the electron screening k. The 
solid gray line is derived from MD simulations [35-38]. The 
dotted black line is from a recent MD simulation of an ultra- 
cold plasma [10]. The results of this work are also plotted as 
black circles with estimated error bars. The two rightmost 
experimental data points correspond to plasmas with low ini- 
tial electron temperatures, as described in the text. Under 
these conditions, the model of Eq (9) tends to overestimate 
the electron temperature, as suggested by the arrows. 



excellent agreement between this result and the predic- 
tions based on Eq. (19) is readily apparent. 

In Fig. 3 we also plot our experimental results. The ex- 
perimental determination of the electron and ion temper- 
atures is described in Sec. HI. There is excellent agree- 
ment between the experimental data and the two simula- 
tions described above. The rightmost experimental data 
point is measured in a plasma evolving from a Rydberg 
gas excited « 10 cm~^ below the ionization potential. 
The second rightmost point corresponds to a plasma ex- 
cited right at threshold. The expansion model we used 
to find the electron and ion temperatures for these two 
plasmas tends to overestimate the electron temperature^ 
at the very early times [34]. Thus we would expect the 
actual electron temperature to be; lower, corresponding 
to larger values of 5 clS suggested by the arrows in the 
plot. 

The data in Fig. 3 demonstrate the validity of the as- 
sumptions used in deriving Eq. (19). The data show that 
electron screening substantially reduces the ion tempera- 
ture, resulting in increased values of F^j. Electron screen- 
ing significantly mitigates the effects of DIH, which is the 
source of ion heating at these early times. Our lowest 
temperature plasmas have F^^ = 4. 
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C. Screened potential energy vs. screened ion 
temperature 

In non-ncutral plasmas, the parameter F^^ completely 
defines the ion-ion interactions. However, in neutral plas- 
mas, a ion-ion interaction necessarily includes contribu- 
tions from the electrons. When the shielding length be- 
comes comparable to the distance between ions, when 
K — >■ 1, the relevance of F^j is questionable. 

One might be tempted to look at Eq. (3) and assume 
that the "effective" coupling constant is F = F^^ exp (— k). 
However, that would overestimate the influence of screen- 
ing. For small k, corresponding to the limit of weak 
screening, the first-order correction in that model would 
be linear in k. This is clearly not the case, as MD sim- 
ulations show. For example, the ion temperature and 
density at the liquid-solid phase transition clearly has no 
linear term [see Fig. 1 in Ref. 37]. 

The idea of calculating F^^ is somewhat problematic in 
neutral plasmas. The F parameter is supposed to repre- 
sent the ion-ion nearest neighbor potential energy divided 
by the ion temperature. The problem arises because the 
ions and electrons are also correlated, and F^^ becomes 
important [see Eq. (10)]. There is a potential energy as- 
sociated with that is shared by both the electrons and 
the ions. Similarly, because the electrons follow the ions, 
there is also a F^^ term that becomes important and that 
will mimic the F^^ behavior. When trying to calculate the 
screened ion coupling parameter, F, it is not immediately 
clear which potential energy is appropriate to include in 
the calculation. They are all important and they all are 
connected to the ion density and temperature. 

This distinction is important to make. The thermody- 
namic properties of non- neutral plasmas depend on F-^. 
These properties can be translated into the realm of neu- 
tral plasmas with the idea that weak electron screening 
modifies them only slightly. However in ultracold neu- 
tral plasmas where k = 1 is achievable, the F,-; scaling of 
these properties is not immediately clear. This is partic- 
ularly the case when the F-^ is determined by k, such as 
we show in Fig. 2. 

In light of the fact that all of the electron and ion cou- 
pling parameters arc important and interconnected, we 
can simply define F to be the total potential energy of the 
system divided by the kinetic energy of the ions. Tak- 
ing U as the total potential energy and K to be kinetic 
energy of the ions, we write 

f = ^ = ^ , (20) 

where the conservation of energy is, trivially, Uj + Kj = 
Up + Kp and we have assumed that the electron tem- 
perature docs not change from the initial to final state. 
Because the ions start out with mK temperatures, we 
can set Kj = Q. Summing up the contributions of the 
electrons and ions to the total initial potential energy 
gives Uj = Q. This can be seen in two ways. One is that 



the initial state is completely uncorrelated and neutral 
and therefore the total potential energy must be zero. 
The other is to argue that the electron-ion potential en- 
ergy terms are negative and exactly cancel the electron- 
electron and ion-ion potential energy terms. Either way, 
we end up at the conclusion that the magnitude of the 
screened coupling parameter is 

f = l. (21) 

Even though electron screening reduces the ion temper- 
ature (see Fig. 2), it reduces the potential energy by 
exactly the same amount so that the ratio of potential 
energy to kinetic energy is always 1. 

This result, Eq. (21), is true for all neutral systems 
in which there is no external source of heat for the elec- 
trons and when there is no correlation in the initial state. 
The agreement between the experimental data and MD 
simulations in Fig. 2 suggest that three-body recombi- 
nation and electron- Rydberg scattering have not signifi- 
cantly increased the electron temperature at these early 
times, because those heating terms are not included in 
the MD simulation. If the electrons are heated, then the 
potential energy Up in Eq. (20) goes down and F will 
increase. 

We note that the final state of the plasma cannot be 
completely determined by energy conservation alone be- 
cause of the two-temperature nature of the UNP. For a 
given initial energy, there are many possible values of the 
final temperatures and that correspond to a correct 
final energy, at least in principle. Of course, if a true equi- 
librium state could be reached, the plasma would have 
Tj = Tg and the final state would be deterministic. This 
suggests that more work on the quasithermodynamics of 
two-temperature plasmas is warranted [40, 41]. 

VI. CONCLUSION 

In this paper we present experimental measurements of 
laser-induced fluorescence from the ions in an ultracold 
neutral plasma. Fluorescence measurements are made 
when the probe laser frequency is scanned over the emis- 
sion lineshape. From these fluorescence data we extract 
the rms velocity distribution as a function of time by fit- 
ting the data to a Voigt profile. An expansion model is 
used to find the electron and ion temperatures. This is 
done over a range of initial electron temperatures, which 
allows us to study the effect of electron shielding on ion 
equilibration at early times. Information about the ion 
and electron temperatures is used to calculate the strong 
coupling parameter F^j and the electron shielding k. We 
compare our experimental results with molecular dynam- 
ics simulations and theoretical calculations for the ion 
strong coupling in ultracold plasmas as a function of the 
electron shielding. We find that our experimental data 
show good agreement with MD results. We generate plas- 
mas with very cold ions because electron screening miti- 
gates the effects of DIH. However we also find that elec- 
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tron shielding softens the ion interaction strength, which 
has the net effect of keeping the ratio of potential energy 
to kinetic energy constant for all values of k. 

Our results indicate that it may be possible to use elec- 
tron screening to generate a strongly coupled plasma with 
> 4. This could be done by ionizing a low-density 
atom cloud with very low initial electron temperature. 
The low density will reduce the time scale for electron 
heating due to three-body recombination. The plasma 
electrons could be heated so that the ions are adiabat- 
ically shifted into their equilibrium positions in an un- 



screened plasma as n is reduced to zero. A large initial 
size for the plasma would also reduce the time scale for 
the plasma to expand. 
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